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We explore charged-Higgs CP-violating effects in an intriguing two-Higgs 



5_i ■ doublet model which accords special status to the top quark. In this model 

the heaviness of the top quark originates naturally from the much larger VEV 
of the second Higgs doublet compared to that of the first. The phenomenology 
of this model is quite distinct from that of the usual formulations of the two- 
Higgs doublet model. In particular, the model can easily account for the 
observed CP violation in the kaon sector even if the CKM matrix is real. The 
associated non-standard CP phase can be monitored through measurements 
of the time-dependent CP asymmetry in B — > ipK$ in experiments at the 
upcoming i3-factories. 
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Despite the stunning successes of the Standard Model (SM) it is widely believed that it 
cannot be the complete theory. It is quite possible that the extraordinary mass scale of the 
top quark is giving us a hint as to the nature of the physics beyond the SM. In this note we 
pursue this theme within the context of an atypical two-Higgs doublet model (2HDM) which 
gives the top quark a special status. The gigantic mass scale of the top quark arises naturally 
in this extension of the SM since the top is the only fermion receiving its mass from the 
vacuum expectation value (VEV) of the second Higgs doublet. All the other fermions receive 
their masses (mostly) from the VEV of the first Higgs doublet. Besides providing a natural 
explanation for the largeness of mt, the introduction of the second Higgs doublet has an 
important bonus in that it allows a new source of CP violation: the charged Higgs sector of 
the model contains a CP-violating phase ( u 8" ) in addition to the usual Cabibbo-Kobayashi- 
Maskawa (CKM) phase |J of the SM. This phase may play a useful role in baryogenesis 
since the CKM phase is believed not to be able to account for the baryon asymmetry of the 
universe in such calculations |2] . Perhaps the most exciting feature of the new phase S is that 
it can be monitored in experiments at the upcoming B factories through time-dependent 
CP studies of B -> ^K s . 

The specific model which we consider was introduced in Ref. || and is essentially a special 
case of the 2HDM of type III (i.e., Model III) Q. We regard this model as an effective low- 
energy theory in which only the (slightly extended) Higgs sector and the regular SM particles 
remain as relevant dynamical particles. The above-mentioned ansatz regarding the quark 
masses leads to a very specific pattern of flavour- and (in general) CP-violating couplings. 
In particular, this model exhibits tree-level flavour-changing neutral Higgs interactions, but 
only in the up-type quark sector. It was shown in Ref. |3j that this model can have significant 
effects on the electron electric dipole moment and on D-D mixing. 

The main purpose of this work is to point out that this top-two-Higgs-doublet model 
(T2HDM) can give rise to a very significant amount of non-standard CP violation in charged 
Higgs interactions ||. Contrary to conventional wisdom, perhaps, these non-standard inter- 
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actions are not particularly well-constrained by flavour- changing neutral current processes. 
For the physics of interest to this work the important unknown parameters of the model are 
tan/5 (i.e., v 2 /vi), m R and 5. We have delineated the allowed parameter space of the model 
by considering the repercussions for K and B decays and for K-K and B-B mixing. Since 
the parameter space is relatively large, we have chosen to focus on three "case studies." In 
the first we take the SM CKM phase 7 || to be identically zero. This scenario is particularly 
interesting since the CP-violating phase 6 in the Higgs sector of the theory is then entirely 
responsible for the CP violation observed in the kaon system. In the second case study 
we set 7 equal to 68°, which is the central value of the SM fits (i.e., 7 = 68° ± 15° |7j). 
For the third case, 7 is taken to be negative and equal to —45°. The important point is 
that, in general, the CP asymmetry for B — > ipKs in this model is appreciably different 
from the value expected within the context of the SM. One feature which distinguishes the 
present analysis from many previous studies of non-standard effects in B — > ipKs is that in 
the present case the effect comes from the new CP-violating phase in the B and B decay 
amplitudes, not from a new phase in B-B mixing || . 

Let us briefly recapitulate some important features of the model of Ref. ||. Consider 
the following Yukawa Lagrangian: 

C Y = -L L <h.El R - Q L <t>xFd R - Q^GI^ur - Q L ^ 2 Gl^u R + h.c, (1) 

where 0« = i<7 2 4>* (i = 1,2), and where E, F and G are 3x3 matrices in generation space; 

= diag(l, 1, 0); 1^ = diag(0, 0, 1); and Ql and Ll are the usual left-handed quark and 
lepton doublets. This Lagrangian gives special status to the top quark, as evidenced by 
the fact that only 02 (and not </>i) couples to t R . Let us set the VEVs of <p\ and <p2 to be 
V\/y2 and V2e 1<T / 'a/2, respectively. In keeping with the spirit of this model, we require that 
tan/3 = v 2 /vi be relatively large (say at least of order twenty). 

The expression for the quark- charged- Higgs Lagrangian takes the form: 

= ( g /y/2m w ) {-u L Vm d d R G + - tan(3H + +u R m u Vd L G + - tan j3H + 

+u R T)Vd L [tan + cot 0\ H + + h.c. } , (2) 
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where G ± and H ± represent the would-be Goldstone bosons and the physical charged Higgs 
bosons, respectively. Here m u and are the diagonal up- and down-type mass matrices, V 
is the usual CKM matrix and £ = m u U R \^U r. U r is the unitary matrix which diagonalizes 
the right-handed up-type quarks. Since £ is in general not a diagonal matrix, the last term 
in Eq. (Q) can give rise to unusual couplings that can violate CP in non-standard ways. To 
see this, let us write the unitary matrix Ur in a form similar to that used in Ref . || : 

\ 



U 



R 




(3) 

Inserting Ur into 



where e ct = m c /m t and where £ is a complex number of order unity 
the definition of E, we obtain 

/0 \ 

S= m c e^|e| 2 m c e ci r y/l - |e rf £| 2 . (4) 
V m c ^l - \e c £\ 2 rrn (1 - |e ct £| 2 ) J 

It is clear from this expression that this model generically has a rather large non-standard b-c- 
H vertex. This can lead to significant CP violation effects in B — > ipK$- Consider, in partic- 
ular, the time-dependent CP-asymmetry a{t) = \r(B(t)) - T(B(t))] / \r(B(t)) + T(B(t))], 
for B — > ipKs- In the SM a(t) is free of hadronic uncertainties and is given by osm(^) — 
— sin (2/3ckm) sin (AMt) |§, where AM = Mr h — Mr l is the mass difference between the 
neutral B mesons and /3qkm = arg ( — KdKb/^^tb)- Existing experimental information 
seems to constrain /3ckm quite tightly: sin(2/3cKM) = -75 ± .10 0. The asymmetry a(t) is a 
very clean way to measure the CKM angle /?ckm and thereby provides an important test of 
the SM [0 



CP-odd effects due to new physics can affect the asymmetry a(t) in several qualitatively 
different ways. The most direct effect comes about if we allow the new interactions to be 
partly (or fully) responsible for the observed CP violation in the kaon system. In this case 
Pckm itself can be different from the value expected within the context of the SM, and can 
even be zero or negative. For a given value of Pckm there can also be new contributions 
to a(t) arising from B-B mixing and from CP violation in the B and B decay amplitudes 



themselves. The B-B mixing effect, being suppressed either by 1/ tan 4 /? for large tan/3 or 
by ml/m 2 H , is generically quite small in this model and may safely be ignored. By way of 
contrast, CP violation in the amplitudes can in principle be appreciable in this model, leading 
to a sizable correction to the asymmetry a(t) compared to the value expected within the 
SM. Ignoring perturbative QCD effects, there are two tree-level diagrams which contribute 
to B — > ipKs: the SM PF-mediated diagram and the new charged-Higgs-mediated diagram. 
The effective Lagrangian is then well-approximated by 

C cS ~ -2V2G F V cb V c * s [c L ^b L s L Yc L + 2(e i5 c R b L s L c R ] + h.c. , (5) 

where (e lS = (l/2)(l^ b /V^)(m c tan/3/m#-) 2 £* with ( taken to be real and positive. In writing 
the above expression, we have neglected terms which are subdominant for tan/3 3> 1 and 
|£| = 0(1). Note that the Higgs-mediated contribution can in principle be large since it is 



proportional to V^/V^,~25 as well as to tan (3 |TT |. 

Using Fierz identities and assuming factorization, we may write the total amplitude for 
B — > ipK$ as A = A(B — > ipK$) — Asm 1 — C e ~ lS • To the extent that factorization holds, 
there is no relative strong phase between the W- and charged-Higgs-mediated diagrams, so 
that the B and B decay amplitudes have the same magnitude and differ only by a CP- 
violating phase. We may thus write A/ A = (Asm/ Asm) exp(— 2i^), where d represents the 
correction to /3ckm due to the Higgs-mediated contribution to B — > ipK s ; i.e., the time- 
dependent asymmetry a(t) now measures /?ckm + $ instead of Pckm- The angle i? is simply 
given by tan$ = ( sin 6/(1 — (cos 5) and can, in principle, be rather large. 

Let us discuss some of the experimental constraints on The most stringent experimen- 
tal constraints on the parameter space of the T2HDM come from b — > 57 and K-K mixing. 
b — > 57 places an important constraint on the parameter space of the T2HDM just as it 
does for Model II, where a tan /3-independent lower bound of approximately 370 GeV may 
be placed on the charged Higgs mass ||12|| . The situation in the present case is somewhat 
complicated by the fact that the Higgs contribution in this model has a CP- violating phase. 
This extra degree of freedom means that the new contribution to b — ► 37 is complex and 



need not interfere constructively with the SM contribution at LO, as is the case in Model II. 
As a result, the shape of the excluded region in the tan/?-m# plane is strongly dependent 
on the CP-violating phase S. 

In order to minimize the theoretical uncertainties associated with b — > 57, it is customary 
to consider the ratio R = Bib — > X s j)/B(b — > X c eV). The branching ratio appearing in the 
numerator of this expression has been investigated by both the CLEO II and ALEPH col- 
laborations . The experimental results may be combined to obtain the weighted average 
R ex P = (2.90 ± 0.46) x 10~ 3 . Considerable progress has also been made on the theoretical side 



of this decay in recent years. In particular, complete NLO QCD corrections JT4| as well as 
some two-loop electroweak corrections [O] have now been incorporated into the calculation 



of R within the context of the SM 16 



The T2HDM can affect the ratio R through diagrams similar to the usual SM dia- 



grams [17] in which the W boson is replaced by a charged Higgs. The Higgs contribution 
gives a correction to the matching condition for the LO Wilson coefficient C^p at the scale 
my/ so that the relevant effective Wilson coefficient at the scale fi ~ 0(m,b) is modified 
according to C^ cS (fi) — > C^ eS (fi) + (a s (m w ) / 'a s (/i)) 16//23 '5C?(mw) , where 



5Cf\m w )=Y J K u [~ tan 2 (3 + i^ T V*) us (tan 2 + l) /m u V* s ] (6) 

u=c,t 

x {B(y u ) + A{y u ) [-1 + (sV)^ (cot 2 /? + l) /m u V ub ] /6} . 
In this expression k u = ±1 for u = c,t; y u = (m u /mH) 2 ', and A and B are standard 



expressions [ 18|Jl^| . 

In order to compare our theoretical calculation of R theor y with the experimental number 
quoted above, we assign a theoretical uncertainty of 8R theoiy = 0.7 x 10~ 3 to R theovy M 



This uncertainty corresponds roughly to the variation of R thcor y as fi is varied between m b j2 
and 2171b. Adding the experimental and theoretical uncertainties in quadrature, we find the 
following constraint on R thc ™y : 0.0012 < R thcory < 0.0046. 

We next consider the implications of the T2HDM for K-K mixing. The total short 
distance contribution to Arrix is given by 
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(Am K ) SD = (G 2 F /6n 2 )f 2 K B K m K \ 2 c x [m 2 cVl + (m 4 tan 4 P/Arr^M] (7) 

where the first term is the usual SM contribution and the second term is the dominant 
contribution in the T2HDM due to the HHcc box diagram. Here Bk is the usual bag factor, 
A c = VcsV* d , and r)i and r)[ are the QCD corrections to the two box diagrams. The SM top 
quark contribution is a few percent of the charm quark contribution and is not included 
here. Similarly, the contributions from WHcc and other box diagrams are negligible in the 
large tan f3 limit and are ignored. 

In order to numerically deduce the allowed parameter space subject to the Attik con- 
straint we use the method described in J7|. Assuming the magnitude of the long distance 
contribution to Am^ to be no larger than 30%, we find the 95% CL. limit m#/ tan 2 (3 > 0.48 



GeV for tan/3 > 10 [20|. Note that the tan 4 /? enhancement in Am^ generally leads to a 
very severe lower bound on the Higgs mass for large tan/? - this is a unique feature of the 
T2HDM. 

Due to the tan 4 /? dependence, the CP- violating parameter €k also receives its largest 
correction from the HHcc box diagram in the large tan/? limit. The SM contribution is well 



known pl |. The dominant Higgs contribution is given by 



,i? i^n D A \4_/ 

~-K 



e^C E B K AXXJp 2 + r ] 2 sm^ + 6)\^\(m c t8.nf3)yAm 2 w m 2 H (8) 



where A, A, p, and 77 are the CKM parameters in the Wolfenstein parameterization p2 
7 = tan^ 1 r]/p is the CKM phase @, and C € = G 2 F f 2 K m 2 w m K /Q^ 2 Am K = 3.78 x 10 4 . 

Since 7 is essentially a free parameter in this model ||23|| , we can obtain bounds on the 
parameter Y = sin(7 + 5)|^|(tan^/20) 4 (200 GeV/m H ) 2 for any given value of 7 by allowing 
V P 2 + V 2 t° vary within its la uncertainties derived from b — > ucv. (We also allow A, etc. 
to vary, as above.) For 7 = 0° we obtain the 95% CL bound 0.08 < Y < 0.39; the bound 
becomes 0.14 < Y < 0.65 for 7 = —45°. If we assume that 7 takes its SM central value 
of 68° [0, the bound becomes —0.085 < Y < 0.08. Unlike the constraint on (m#, tan/3) 
coming from Amx, the one coming from ex depends on |£| and S. 
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Figure 1 shows plots in the tan (3-rriH plane for a few representative values of the param- 
eters of this model. In all of these plots we have fixed |£| to be unity, which is a natural 
choice in this model. The shaded bands in Fig. 1 correspond to the regions allowed by all 
three constraints. The amplitude of the time-dependent CP asymmetry in B — > ipKs is now 
given by a^ Ks = sin [2(/3 CK m + #)]■ Contours of constant i? have a very simple form in the 
tan (5-iriH plane: they are lines of constant slope emanating from the origin of that plane. 
We have suppressed such contours to avoid overcrowding the plots. Figure 2 shows the 
allowed values of sin [2(/?ckm + $)] as a function of the CP-violating parameter 5. In order 
to clearly illustrate charged Higgs effects, we fix a/ p 2 + f] 2 = 0.41 (this is its central value 
0) so that /3qkm is uniquely determined for a given 7. This situation might well correspond 
to a future scenario in which measurements of b — > uev have become more precise. 

For illustration let us focus on three representative choices for 7 (see Fig. 2). The first 
corresponds to a real CKM matrix; i.e., 7 = 0°. In this case the phase 6 is solely responsible 
for the observed CP violation in the K-K system. For a given 5, the upper limit of the 
allowed region corresponds to the intersection of the b —>■ 37 curve with the more stringent 
of the Am>K and the ex curves in the tan f3-mn plane. The excluded regions for 6 come 
dominantly from the €k constraint. It is interesting to note that small values (~ 5°) for 
S are not ruled out by 6% due to the tan 4 /? enhancement factor in Eq. @. This scenario 
gives a very clear signature: instead of measuring an asymmetry of order sin2/?cKM ~ 0.75 
as expected in the SM, one would find an asymmetry no larger than about 0.27. 

Figure 2 also shows the allowed range of the CP asymmetry a^K s for another phenomeno- 
logically interesting case, 7 = 68°. It is evident from this plot that significant deviations 
from the pure SM value (i.e. the horizontal line corresponding to $ = 0) are possible in this 
scenario. Depending on the value of 5, the partial rate asymmetry in B — > ipKs can range 
between ~ 42-90% and therefore may be appreciably different from the SM expectation of 
~ 75%. 

As the third scenario, we show in Fig. 2 how the CP asymmetry a^K s can even have 
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an opposite sign relative to the standard model expectation. This generally occurs when 
7 takes negative values which in turn gives a negative Pckm- As an example, the allowed 
region in the a^K s -8 plane is shown for 7 = —45°. This scenario gives a distinct signal for 
physics beyond the SM. 

To summarize, we have presented a study of CP violation in a two-Higgs-doublet model 
that accords a special status to the top quark and accommodates the remarkable heaviness of 
the top quark rather naturally. While the model has repercussions for numerous experiments 
24j, we have focused here on its effects on the time-dependent CP asymmetry in B — > ipKs- 



This asymmetry is of particular interest due to the fact that reliable theoretical predictions 
can be made here and also because intense experimental activity by the B factories is 
anticipated on this decay mode in the very near future. The asymmetry in this model can 
differ significantly from the predictions of the SM and can even change sign with respect 
to the SM expectation. These features may be helpful in leading to the discovery of new 
physics. 

We thank D. Atwood, A. Czarnecki and O. Vives for helpful discussions and comments. 
This research was supported in part by the U.S. Department of Energy under contract num- 
bers DE-AC02-98CH10886 and DE-FG02-91ER40681 (Task B). K.K. was also supported in 
part by the Natural Sciences and Engineering Research Council of Canada and by an SRTP 
grant at Taylor University. G.W. is grateful for the hospitality of the High Energy Theory 
Group at Brookhaven National Lab and of the Fermilab Theoretical Physics Department. 



9 



REFERENCES 



[1] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 
[2] See, e.g., A.G. Cohen et al, Ann. Rev. Nucl. Part. Sci. 43, 27 (1993). 
[3] A. Das and C. Kao, Phys. Lett. B 372, 106 (1996). 

[4] See, e.g., W.-S. Hou, Phys. Lett. B 296, 179 (1992); M. Luke and M.J. Savage, Phys. 
Lett. B 307, 387 (1993); D. Atwood, et al, Phys. Rev. D 55, 3156 (1997); M. Sher, 
|hep-ph/980959q . 

[5] Assuming CP violation arises from the Higgs potential and using approximate global 
family symmetries to suppress FCNC, extra phases can also appear in the charged Higgs 
couplings in a 2HDM (model III). Consequences of this type of model have been studied 



by Y. L. Wu and L. Wolfenstein, Phys. Rev. Lett. 73, 1762 (1994), and Y. L. Wu, [hep 



ph/ 940424 1| . One crucial difference from the T2HDM in which we are interested is that 
tree-level contributions to KK and BB mixings are absent in the T2HDM but are 
generally present in the former. 
[6] See, e.g., the article by H. Quinn in Review of Particle Physics, C. Caso et al, Eur. 
Phys. J. C 3, 1 (1998). 



[7] S. Mele, |hep-ph/98084TT1 ; see also P. Paganini et al, |hep-ph/ 9 7112611 



[8] See, e.g., M. Gronau and D. London, Phys. Rev. D 55, 2845 (1997); Y. Grossman and 
Y. Nir, Phys. Lett. B 313, 126 (1993); Y. Grossman and M. Worah in Phys. Lett. B 
395, 241 (1997). 

[9] The form which we have chosen for Ur differs slightly from that used in Ref. ||, the main 
difference being that we neglect possible mixing between the first and third generations. 
If this mixing is suppressed by m c /m t or more, then it has a negligible effect in the 
processes which we will consider. Dropping this term also allows us to write Ur in a 
manifestly unitary form. 

[10] The subscript is added to /3qkm to distinguish it from f3 = tan _1 (t> 2 /t>i). 

[11] Indeed this non-standard source may contribute ~ 20% to the observed branching ratio 



10 



for B — > which is much smaller than the intrinsic theoretical uncertainties in the 
absolute rate calculation. 
[12] See, e.g., M. Ciuchini, et al, Nucl. Phys. B527, 21 (1998). 

[13] M.S. Alam et al. (CLEO Collab.), contribution to ICHEP conference, Vancouver, 
Canada, July 23-29 1998; R. Barate et al. (ALEPH Collab.), Phys. Lett. B 429, 169 
(1998). 

[14] K. Chetyrkin et al, Phys. Lett. B 400, 206 (1997); erratum ibid B 425, 414 (1998); 

A.J. Buras et al, Phys. Lett. B 414, 157 (1997) and [hep-ph/9707482 . 
[15] A. Czarnecki and W.J. Marciano, Phys. Rev. Lett. 81, 277 (1998). 
[16] Since we restrict ourselves to computing the LO effect of the new Higgs interactions, 

for consistency, we also incorporate only the LO SM contribution to b — > sj. 
[17] T. Inami and C.S. Lim, Prog. Theor. Phys. 65, 297 (1981). 
[18] B. Grinstein et al, Nucl. Phys. B339, 269 (1990). 

[19] Note that we are neglecting here the scalar operator Os = c^jjs^c^, which receives 
contributions at tree-level in this model [see Eq. (|^)]. This term could potentially con- 
tribute to C^ eS (/j,) through operator mixing and our neglect of it is an approximation. 
Analogous terms are also present in Models I and II, but are safely neglected there 
because the Wilson coefficient Cs{mw) is generally very small in these models. 

[20] We take all error distributions (both for the theoretical and the experimental quantities) 
to be rectangular, with widths ±ler. 

[21] G. Buchalla, et al, Rev. Mod. Phys. 68, 1125 (1996). 

[22] L. Wolfenstein, Phys. Rev. Lett. 51, 1841 (1983). 

[23] Note that the choice 7 = is disfavored at the la level by the combined constraints 
from B-B mixing and b — > uev |7j. However, to consider the full range of possibilities 
for this model - and given the theoretical uncertainties implicit in these constraints - 
we choose to adopt a rather conservative approach with respect to their interpretations. 

[24] For example, the charged Higgs exchange contribution to the neutron EDM in this 



11 



model can be as large as 10 27 — 10 26 e cm. 



12 



FIGURES 
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FIG. 1. Experimental constraints on the T2HDM from b — ► 57, e^; and Am^. The allowed 
regions are shaded. The dots indicate where the largest values of occur: these are 1? = 7.0° and 
•d = 5.5° in (a) and (b), respectively. The plot for 7 = —45° is similar to that of 7 = 0°, and is not 
shown. 
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FIG. 2. Amplitude of the time-dependent CP asymmetry cl^k s = sin(2/3cKM + 20) versus the 
non-standard CP-odd phase 8 for B — > ^ifs. The top horizontal line is for the SM assuming the 
best fit value (sin2/?cKM = 0.75) of Ref. [6]. The shaded regions correspond to the allowed ranges 
of the asymmetry in the T2HDM for three representative choices of 7: 7 = 68° is the best fit of 
Ref. [6], 7 = 0° corresponds to a real CKM matrix, and 7 = —45° changes the sign of ol^k s relative 
to the SM expectation. 



14 



